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The spectral flow in the cores of moving dislocations is found to play an important role in the 
dynamics and transport of charge density wave (CDW) by significantly modifying the hydrodynamic 
effective action of the condensate. The analogy of spectral flow in CDW with the baryogenesis in 
the early universe is pointed out. 
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The dynamics of the dislocation lines in charge density 
waves (CDW's) jlj is one of the key processes which gov- 
ern the sliding (nonlinear) conductance of the incommen- 
surate CDW with three dimensional (3D) order. When 
current contacts are attached to the side of the chains, 
a sliding CDW implies that wave fronts are added or re- 
moved at the contacts by the nucleation or annihilation 
of dislocation loops as shown in Fig. [I]. The narrow 
band noise @J3] and the nonlinear current- voltage char- 
acteristics Pppf have been attributed to this phenomena. 
Although this idea explains most qualitative features of 
CDW very well, the dynamics of the dislocations is not 
yet well understood. Especially, the gauge invariance or 
charge conservation condition of CDW has not been in- 
vestigated seriously. 

In this Letter, we first show that the gauge invariance 
of the conventional hydrodynamic effective action of the 
condensate is violated in the presence of moving disloca- 
tions. The origin of this is shown to be the neglect of the 
"spectral flow" in the cores of the moving dislocations. 
By explicitly taking the effects of the spectral flow into 
account, we derive a gauge invariant action which cor- 
rectly describes the dislocation dynamics. We find two 
correction terms arising from the spectral flow, namely, 
the correction to the interaction of the condensate with 
the electromagnetic field and to the condensation energy. 
The nucleation process of the dislocation loop is analyzed 
based on the obtained action and the nonlinear conduc- 
tion of the CDW is discussed in the light of our findings. 
Finally we discuss a remarkable similarity of the dislo- 
cations in CDW with the so-called cosmic strings in the 
early universe, suggesting the intriguing possibility of the 
CDW system as a theoretical and experimental model for 
cosmological field theories. 

First we discuss the gauge invariance of the effective ac- 
tion in the presence of moving dislocations. It has been 
widely accepted that the condensate can be described by 
a hydrodynamic action, which is valid on a length scale 
larger than £ Q 's B, with £ a being the amplitude coher- 
ence length (a = x,y, z depending on the direction). The 
action consists of the elastic part S c \ and the interaction 



part with the electromagnetic field Smt, 

S — S e \ + Sint 

drdr- 



- {c - 2 (d T 9f + (d x 0) 2 
S^=J dTd viL{ V d x e + AAB) 
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where 8(t, r), ip(r, r) and A x (t, r) are the phase of the or- 
der parameter, the scalar potential and the x-component 
of the vector potential, respectively, r = (x,y,z) and r 
is the imaginary time. The conducting chains are paral- 
lel to the x-axis. K — N±f c hvp /2tt and L — eN±f c /ir, 
where f c , Vy, N± are the condensate fraction, the Fermi 
velocity parallel to the chains and the areal density of 
the chains, respectively, cq — i>f y^m/m* is the phason 
velocity with m* /m being the effective mass ratio of the 
condensate to the normal electrons. J y (= £,y/£,x) and 
7z(= ^z/ix) are the anisotropy parameters. 

According to the action, Eq. (||), the charge density 
Pc(t, r) and current density j c (r, r) of the condensate are 
expressed by, 



p c = Ld x e, j c = —L e x d T 9, 



(4) 



where e x is the unit vector parallel to x-axis. The charge 
conservation law then reads, 



l (d T d x e - d x d T e) = r, 



qp- 



(5) 



Although r qp (r, r) vanishes when 8(t, r) is a single valued 
smooth function of r and r, it does not vanish in the 
presence of moving dislocations [M and gives [O 



T qp = -2ttL^ J dl(diR v x d T R u ) x S^(r-K u ), 



(0) 



where S^ir) and R v (t,1) are the d-dimcnsional 6- 
function and the positions of the v-th dislocation line at 
imaginary time r, respectively. I is a parameter along the 
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line. In this Letter we consider only the dislocations with 
2ir phase singularity and those with higher (Att, 6tt • • •) 
singularities are disregarded because of the higher exci- 
tation energies. The Eq. (^) clearly shows that charge 
conservation is violated in the cores of the dislocations, 
and implies the conversion of the quasiparticles into con- 
densate electrons. Although this phenomenon has been 
discussed earlier by several authors PJPH, its effects on 
the effective action of the dislocations have not been con- 
sidered. 

The origin of r qp becomes clearer by examining the 
energy spectrum of the quasiparticles in the core re- 
gion. We consider chains of a finite length l x (0 < 
x < l x = 1) and assume that a straight dislocation 
line parallel to the z-axis is moving from y = to 
y = l y along the plane x — l x /2. (See Fig. || (a).) 
We calculate the energy spectrum of the chain located 
at y = ly/2. The complex order parameter, A(r, r), on 
the chain is parameterized as A(x,Y v ), where Y v is the 
y-coordinate of the dislocation line. The anisotropy of 
the system and the interchain electron hopping is ne- 
glected for simplicity. We take A(x,Y v ) to be constant 
at x = 0, 1, by assuming perfect pinning at the ends 
of the chains. A(x,Y v ) must be quasiperiodic in Y v , 
A(x,Y v ) oc exp(— 2wix) A(x,Y v + l y )), corresponding to 
the addition of a wave front of CDW due to the disloca- 
tion By numerically diagonalizing the Hamiltonian, 

n J ax* y A*(x, Y v ) hv F td x J W ' [t) 

where * = (lZ(x) 7 C{x)) with TZ(x) and C(x) being the 
field operators for right and left moving electrons omit- 
ting the spin degree of freedom, we obtain the energy 
spectrum of the quasiparticles shown in Fig. || (b). It 
clearly shows the spectral flow phenomenon; one energy 
level is transferred from above to below the energy gap 
through the core of the dislocation. When the energy 
level is occupied (vacant) at Y v = 0, a quasielectron in 
the conduction band (quasihole in the valence band) is 
deleted (created) during this process. Since this process 
is strongly localized to the cores of the dislocations, the 
size of which is given by the amplitude coherence length 
£ Q 's, it cannot be described by the conventional hydro- 
dynamic effective action of the condensate which is valid 
on a scale larger than £ a 's. 

Our next task is to derive the effective action of the 
condensate which correctly describe the spectral flow in 
the cores of dislocations. We found that two correction 
terms to Eq. (|^) are needed. One is the correction term 
to the electromagnetic interaction, Si n t, and the other is 
related to the change of the total condensation energy 
by the spectral flow. The former should be introduced 
to compensate the violation of the gauge invariance dis- 
cussed above. We first focus on a single chain to ana- 
lyze this situation. The motion of a dislocation is then a 



phase slip process localized at the core |12j]. The quasi- 
particles are created at this position and diffuse over the 
entire chain. The number of created quasiparticles is 
two (including the spin degree of freedom) at K, and 
decreases at higher temperatures in proportion to the 
equilibrium condensate fraction / c . The actual diffusion 
process of the quasiparticles is governed by the mobility, 
which depends on the material. In a dirty CDW sample, 
the motion of the quasiparticles is diffusive, whereas in 
a clean one it can be ballistic. The situations are clas- 
sified by the relation between the quasiparticle velocity 
v qp and the phason velocity cq. The latter characterizes 
the velocity of the propagation of deformations in the 
condensate. The following two limiting cases are useful. 

Clean limit (v qp 3> Co) — Since the effective mass of 
the quasiparticles near the energy gap can be approxi- 
mated by m c g ~ Ao/wp, where Ao is the equilibrium 
energy gap, an applied field accelerates the quasiparti- 
cles to v qp = y eAo/™ ff ~ vf without impurity scatter- 
ing. (The particle- hole process limits the speed.) Since 
vf 3> co the quasiparticles move much faster than the 
condensate in this limit. Then the spatial change in the 
quasiparticle current density at the dislocation position 
can compensate r qp in Eq. (^|) and the correction can 
be implemented by introducing the following term to the 
effective action, 

SSint = -iL J drdr A x (t, r) J dx' Y qp (r, x' , y, z) . 

(8) 

Dirty limit (v qp <C Co) - With impurity scattering, 
v qp is given by —^ qp E x , where E x and /i qp are the applied 
field and quasiparticle mobility, respectively. When v qp 
is much smaller than Co it is a good approximation to 
neglect the motion of the created quasiparticles. Then 
the correction term becomes, 

6S int = %L J drdr <p (r, r) / dr' T qp (r', r) . (9) 

For example, if the mobility is of the order of ~ 
lem 2 /V • s, which is a typical value for the conventional 
samples, i> qp is smaller than co unless E x > 10 6 V/cm, 
which is much larger than the experimental value. The 
dirty limit therefore applies to most of conventional CDW 
samples. 

The second correction term arises from the conden- 
sation energy. The dislocation motion in CDW al- 
ways produces quasiparticles from the condensate. If 
some quasiparticles already exist near the dislocation 
core, the dislocation can convert some to the conden- 
sate and gain condensation energy. Employing Eq. (|^), 
the contribution of the condensation energy to the ef- 
fective action of the dislocation is estimated as SS c d ~ 
± J drdr(Ao/e) J dr'r qp (T', r), where ± depends on the 
situations described above. In most CDW materials, 
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Ao/e is of the order of 0.1 V and is important, espe- 
cially when the x-component of the applied electric field 
is small (for example, in a transverse field condition [H). 

Next we discuss the effective action of the disloca- 
tions taking into account the correction derived above. 
It is convenient to transform here to an isotropic co- 
ordinate system by x = (xq, x\, x%, X3) with xq = cqt, 
xi = x, x 2 = y/j y , x 3 = z/~f z , K = Kjyjz/co and 
L = Ljyjz/co. We express 8^9 in terms of the dis- 
location coordinate since d^9 is well-defined and sin- 
gle valued, although 9 itself is not. As in Eq. (||), 
we parameterize the position of the v-th. dislocation in 
space and time by I and r, respectively. The posi- 
tion of the line element in 4D space is expressed by 
7] v (t,1) = (cot,R„) = {cqt,X v ,Y v ,Z v ). Introducing a 
flux tensor, J- a p(x) [flB) , 



T a0 {x) = 2*J2 dldT (Vu) a {v'u)0^(x - r) u ), (10) 
where r\ v = d T r\ v and rf u = dir\ v , we obtain 

d M 6(x)= [ d^x'e^daG^ix-x'^ix'), (11) 



where G^(x) = -(An 2 
of —dftdft and e^ Qj a 7 is the fully antisymmetric tensor. 
Here the single valued part of 9 is omitted. Employing 
Eq. ( pd| ) we can express all the terms in terms of Rj, (r, I). 
For example, S e \ of Eq. (|3|) is rewritten as, 

S*d = Jd 4 xd 4 x' y G {A \x - x')F aP (x)? af3 (x'), (12) 

where the conservation law of the topological charge of 
the dislocations, d a T a (i — has been used. 

Now we consider the sliding of CDW by the nucleation 
of dislocations near the contacts HHQ. The effective 
action discussed above depends on the electric field in the 
sample, which is conventionally assumed stationary and 
parallel to the chains. The "two fluid approximation" 
yields 



{-d 2 x + \^)E X 



inD- 1 {j qp } x . 



(13) 



where Aq 2 = &e 2 f c / (Uvf), D and {jqp}^, are the 
Thomas-Fermi screening length of the normal state, the 
diffusion constant and the cc-component of the quasipar- 
ticle current, respectively. From this equation we can 
see that the electric field distribution depends on j qp and 
D, in other words, on the density and mobility of the 
quasiparticles. 

The thermal and quantum nucleation rate of a 
dislocation loop in the presence of an applied elec- 
tric field can now be estimated based on the effec- 
tive action derived above. We employ the method 
by Langer and Fisher |l6],[|,[7| and by Duan Q, 



to estimate the thermal and quantum nucleation 
rate, respectively. In the latter case we take the 
bounce solution of a dislocation loop as T)(t, 1?) = 
(cor, X, ^R 2 - (c r) 2 costf, yjR 2 - (c r) 2 sintf), (0 < 
•d < 2tt), and use the instanton method. 

According to the above arguments we can distinguish 
between the clean and dirty samples and between the 
metallic and semiconducting materials, which leads to 
the following picture. 

(1) In clean systems, in which the quasiparticlc mo- 
tion is ballistic, 5Si nt of Eq. (||) must be employed. In 
this case the dislocations behave like an array of elec- 
tric dipoles with moments aligned perpendicular to the 
chains H] which do not couple with E x . Therefore the 
nucleation of the dislocation loop does not lead to a gain 
of electrostatic energy. In addition, since the deviation 
of the quasiparticle density from the equilibrium is small, 
the gain of condensation energy is also negligible. The 
nucleation of the dislocation loops is then unlikely in 
these systems and sliding CDW transport is suppressed. 

(2) In dirty metallic systems, S- m t in Eq. (||) is modified 
by Eq. (|) to 



■Jin 



Si, 



= iL 



-SS h 
d A x 



ip{x) / dx' did' 9(x' ,: 



+ c A x { X )d a 9(x) 



l -^Z„Jdx Jdl (R„ x d{R. u ) x V(X V ), 



(14) 



(15) 



where d\ = d/dxi, d' a = d/dx' and x = (xi, X2, X3). 
In the last equation an approximation cq — > 00 is ap- 
plied, which yields G^'(x) — > — (47r) _1 (5(a;o)|x| _1 and 
the gauge is fixed as 4 = 0. V(X V ) is given by 
J X dxi sgn(xi — X v )E x (x\). Note from Eq. ( |l5| ) that 
S' int is proportional to the area surrounded by the loop. 
In this case E x is approximately constant (~ V/l x ) in the 
entire sample, where V is the applied voltage. 

Since the driving force of the dislocation is strongest 
near x± = and X\ — l x , the nucleation of disloca- 
tion loop is most likely to occur at the edges. The ef- 
fect of the boundary at x± = (or l x ) can be taken 
into account by considering mirror dislocations with re- 
spect to the boundaries, which replaces V{X V ) with 
V b {X v ) = cfei E x { Xl ) (= 2/ ^---) in case of 

X v ~ (~ l x "). The nucleation rate is then found to 
be proportional to exp{ — (Vq/V) 2 } for quantum nucle- 
ation, where Vq = (ir 2 /12)y/m*/m?i 3 "/ y j z Vpf c Nj_e^ 2 , 
and exp(— Vr/V), for thermal nucleation, where Vt = 
irjyJzf c N ± (hv F ) 2 /(ek B T). V rj \V b {X„)\ is the applied 
voltage. 

(3) In dirty semiconducting systems, the same action 
as Eq. ( |l5| ) applies. However, since j qp ~ in Eq. (|13), 
E x in the bulk can be much smaller than V/l x due to 
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the screening by the condensate, which reduces the nu- 
cleation rate. This qualitatively agrees with experiments 
Since the quasiparticle density is also strongly 
modulated, 5S c d becomes important, which causes a 
deviation of the nonlinear I-V characteristics from the 
metallic cases when V < A/e ~ 0.1 V. 

The result obtained in (2) coincides with those ob- 
tained earlier in a different framework, which ap- 
pears to be accidental however. These theories attribute 
the driving force of the dislocations to the elastic stress 
induced in the condensate by the applied electric field, 
which would imply the unlikely scenario that disloca- 
tions nucleate when an electric field is applied without 
contacts. From our theory it is clear that by nucleation 
of dislocations no energy is gained in the above situation, 
since it just results in an over-screening by the conden- 
sate. This demonstrates that our theory describes the 
screening and transport properties correctly. 

Before concluding this Letter, we comment on the sim- 
ilarity of the CDW system with field theories for the early 
universe. A dislocation in CDW is a counterpart of a 
cosmic string, i.e. a topological defect in electroweak 
field; The former creates quasiparticles from the conden- 
sate (Fermi sea) without creating anti-quasiparticles as 
is clear from Fig. ||, whereas the latter produces matter 
without producing antimatter, which is called baryogen- 
csis jjl|[l9| • Although the similarity of baryogenesis with 
the momentum creation in the vortices in superfluid Hc3 
has already been discussed [^0| , the CDW system is more 
similar to the electroweak theory in the sense that it has 
an identical chiral symmetry and chiral anomaly | pl| . We 
therefore expect that further theoretical and experimen- 
tal studies of CDW may shed new light on its field theo- 
retical counterpart in particle physics and cosmology. 

In conclusion, we studied the dynamics of the disloca- 
tions in CDW's taking into account the contribution of 
the spectral flow in the cores. Our theory not only pro- 
vides a natural basis to understand the dynamics of the 
dislocations in CDW's but also indicates new aspects of 
CDW. 
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FIG. 1. A schematic view of the nonlinear conductivity 
measurement. Wave fronts are indicated by shades. Disloca- 
tion loops and their motions are shown. 




FIG. 2. (a) The configuration of the dislocation (©), the 
chains (dotted lines) and wave fronts (shades), (b) The en- 
ergy spectrum of a chain (bold line of (a)) as a function of 
the position of the dislocation, Y v , is shown. All the energy 
levels are doubly degenerated at Y v = 0, l y except for the in- 
nermost two with respect to the energy gap at Y v = 0. In this 
calculation £ = 0.1 is taken. 
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